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ABSTRACT 

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. 

We utilize the bright afterglow of a Population III (Pop III) GRB exploding in a primordial dwarf galaxy as a 
luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first 
heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an 
existing highly resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic 
hydrogen cooling in a halo with virial temperature > I O 4 K. We explore the unusual circumburst environment inside 
the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne 
flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to ( 1 + z) with 
values of about a few cm -3 . In more massive halos, corresponding to the first galaxies, the density may be larger, 
n > 100 cm -3 . The resulting afterglow fluxes are weakly dependent on redshift at a fixed observed time, and may 
be detectable with the James Webb Space Telescope and Very Large Array in the near-IR and radio wavebands, 
respectively, out to redshift z > 20. We predict that the maximum of the afterglow emission shifts from near-IR to 
millimeter bands with peak fluxes from mJy to Jy at different observed times. The metal absorption line signature 
is expected to be detectable in the near future. GRBs are ideal tools for probing the metal enrichment in the early 
IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by 
the first supernova (SN) explosions are likely to reside in low-ionization stages (Cn, Oi, Sin and Fen). We show 
that, if the afterglow can be observed sufficiently early, analysis of the metal lines may distinguish whether the first 
heavy elements were produced in a pair-instability supernova or a core-collapse (Type II) SN, thus constraining the 
initial mass function of the first stars. 

Key words: cosmology: observations - cosmology: theory - galaxies: high-redshift - gamma-ray burst: general - 
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1. INTRODUCTION 

One of the key goals in modern cosmology is to study the 
formation of the first stars and galaxies at the end of the cosmic 

dark ages, and how they shaped the subsequent evolution of the 
universe (Barkana & Loeb 2001 ; Bromm & Larson 2004; Ciardi 

& Ferrara 2005; Bromm et al. 2009; Loeb 2010). The first, 

so-called Population III (Pop III) stars are predicted to have 
formed at - > 20 in minihalos with virial mass M wu ~ 10 6 M 0 

and temperatures T v „ < 10 4 K (Haiman et al. 1996; Tegmark 

et al. 1997 ; Yoshida et al. 2003). The Pop III initial mass function 

(IMF) is thought to be top-heavy (Bromm et al. 1999, 2002; 
Abel et al. 2002), possibly extending to M* > 100 M 0 , but 
recent simulations indicate that the Pop III IMF may be quite 
broad, also including a fraction of lower-mass stars (Stacy et al. 
2010; Clark et al. 2011b; Greif et al. 2011, 2012). The first 
bona-fide galaxies are expected to have formed at a later stage 
in hierarchical structure formation (Bromm & Yoshida 2011), 
when ~ 1 0 8 M 0 halos assembled at z > 10 via the merging 
of progenitor minihalos (Wise & Abel 2007; Greif et al. 2008). 
These systems are often termed “atomic cooling halos,” because 
their virial temperature, T v i r > 10 4 K, exceeds the threshold 
to enable efficient cooling via lines of atomic hydrogen (Oh 
& Haiman 2002). Direct observations of the first galaxies at 
redshifts z > 10 have so far been out of reach. In the coming 
decade, the James Webb Space Telescope ( JWST ) promises to 


directly probe this critical period (Gardner et al. 2006). The 
detection of metal absorption lines in the afterglow spectrum of 
high-redshift gamma-ray bursts (GRBs), imprinted by enriched 
gas in the first galaxies, offers an unusual opportunity to study 
the physical conditions inside them. We may thus be able to 
derive constraints on the temperature, metallicity, ionization 
state, and kinematics in the interstellar medium (ISM) of high- 
redshift galaxies, and in the surrounding intergalactic medium 
(IGM). It is encouraging that such diagnostics can already be 
obtained for bursts at somewhat lower redshifts. An example 
is GRB 081008, where high-resolution spectroscopy with the 
Very Large Telescope has probed the ISM of a host galaxy at 
Z ~ 2 (D’Elia et al. 2011). 

Long-duration GRBs have been shown to be associated with 
the death of massive stars (Stanek et al. 2003; Hjorth et al. 2003; 
Woosley & Bloom 2006). Their high luminosities make them 
detectable out to the edge of the visible universe (Ciardi & Loeb 
2000; Lamb & Reichart 2000; Bromm & Loeb 2002, 2006; Gou 
et al. 2004; Inoue et al. 2007; Meszaros & Rees 2010), with the 
current record held by GRB 090429B at z ~ 9.4(Cucchiaraetal. 
2011). GRBs provide ideal probes of the high-redshift universe, 
including the star formation rate (Totani 1997; Wijers et al. 
1998; Porciani & Madau 2001; Chary et al. 2007; Yiiksel et al. 
2008; Wang & Dai 2009, 2011; Elliott et al. 2012), reionization 
(Gallerani et al. 2008), dark energy (Dai et al. 2004; Wang et al. 
2011), and the IGM metal enrichment (Barkana & Loeb 2004; 
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Totani et al. 2006; Toma et al. 2011; Bromm & Loeb 2012). The 
leading contender for the central engine of long-duration GRBs 
is the collapsar model (Woosley 1993; MacFadyen et al. 2001). 
Because of their predicted high characteristic mass, a significant 
fraction of Pop III stars might end their lives as a black hole, 
potentially leading to a large number of high-redshift GRBs. 
Thus, Pop III stars are viable progenitors of long-duration GRBs, 
triggered by the collapsar mechanism, as long as they can lose 
their outer envelope and retain sufficient angular momentum 
in their center (Bromm & Loeb 2006; Belczynski et al. 2007; 
Komissarov & Barkov 2010; Stacy et al. 2011). It might even 
be possible for Pop III collapsars to occur if the extended outer 
envelope were not lost (Suwa & Ioka 201 1). 

The history of pre-galactic metal enrichment has several 
important consequences for structure formation (Madau et al. 
2001; Karlsson et al. 2012). An early phase of metal injection 
may qualitatively change the character of star formation, from a 
predominantly high-mass (Pop III) mode to a normal, low-mass- 
dominated (Pop I/II) one, once the enrichment has exceeded 
a “critical metallicity” of Z CI j t ~ 1 0 4 Z Q (Bromm et al. 
2001a; Schneider et al. 2002, 2006; Bromm & Loeb 2003; 
Mackey et al. 2003). The transition between these two modes has 
crucial implications, e.g., for the expected redshift distribution 
of GRBs (Bromm & Loeb 2006; Campisi et al. 2011; de Souza 
et al. 2011), for reionization (Cen 2003; Wyithe & Loeb 2003; 
Furlanetto & Loeb 2005), and for the chemical abundance 
patterns of low-metallicity stars (Qian & Wasserburg 2001; 
Frebel et al. 2007, 2009; Tumlinson 2010). It is therefore 
important to explore the topology of early metal enrichment, 
and to determine when particular regions in the universe become 
supercritical (Tornatore et al. 2007; Maio et al. 2010). 

Absorption lines imprinted on the spectra of bright back- 
ground sources, such as GRBs or quasars, are one of the main 
sources of information about the physical and chemical prop- 
erties of high-redshift systems (Oh 2002; Furlanetto & Loeb 
2003; Oppenheimer et al. 2009). These lines are due mainly 
to absorption by neutral hydrogen present in the low column- 
density Ly« forest, and by metals in low-ionization stages (e.g., 
C ii, Si n, Mg ii, Fe ii, O i) which arise in the higher column- 
density gas associated with Damped Lya Absorbers (DLAs). 
The analysis of the spectrum of distant GRB 050904 (Totani 
et al. 2006) has resulted in a wealth of detailed insight into 
the physical conditions within the host galaxy at z — 6.3, and 
Salvaterra et al. (2009) claimed that they identified two absorp- 
tion lines (Si iv and Fe ii), although at poor signal-to-noise ratio, 
in the spectrum of GRB 090423, the most distant spectroscopi- 
cally confirmed burst at z = 8.2 (Salvaterra et al. 2009; Tanvir 
et al. 2009). GRBs as background sources offer a number of 
advantages compared to traditional lighthouses such as quasars 
(Bromm & Loeb 2012). Their number density drops much less 
precipitously than quasars at z > 6 (Fan et al. 2006), and the 
absence of a strong proximity effect, together with the near 
power-law character of their spectra, renders them ideal probes 
of the early IGM. 

In this paper, we discuss the observational signatures of Pop 
III GRBs and study pre-galactic metal enrichment utilizing 
absorption lines in the spectra of high-’ GRBs which were 
imprinted by the first galaxies. Recently, it has become feasible 
to study the formation of the first galaxies, including the metal 
enrichment from Pop III supernovae (SNe), with highly resolved 
cosmological simulations (Wise & Abel 2008; Greif et al. 2010). 
We here place a bright GRB into the simulation box of Greif 
et al. (2010), and derive the spectral signature as the afterglow 


light escapes from the first galaxy, thereby probing the partially 
enriched IGM in its vicinity. The fluxes in the near-IR and radio 
bands may be detectable by the JWST and the Very Large Array 
(VLA) out to z > 20. The structure of this paper is as follows. In 
Section 2, we derive the circumburst density profile of Pop III 
GRBs, followed by a brief description of the underlying first 
galaxy simulation (Section 3). We discuss the properties of 
the GRB afterglow in Section 4, and the metal absorption line 
diagnostics in Section 5, followed by our conclusions. 


2. GRB DENSITY PROFILE 

Two different environments are currently discussed as Pop III 
star formation sites, minihalos and atomic cooling halos 
(Bromm et al. 2009). The latter mode is often termed Pop III. 2 
to indicate physical conditions, such as a higher degree of ion- 
ization, that differs from the canonical minihalo case, possibly 
resulting in somewhat lower masses (McKee & Tan 2008). For 
normal GRBs, either a constant number-density profile, or a 
power-law dependence, as might be expected in the stellar wind 
from the progenitor, is discussed in the literature (Bloom 201 1). 
The GRB afterglow emission sensitively depends on the circum- 
burst density (Ciardi & Loeb 2000; Gou et al. 2004; Inoue et al. 
2007), and this is one of the crucial uncertainties in making pre- 
dictions for Pop III bursts. We next derive the likely circumburst 
conditions in the two possible Pop III host systems. 

2.1. Minihalo Case 


For simplicity, we assume that Pop III GRBs are triggered 
in minihalos close to the cooling threshold for collapse, which 
is only weekly dependent on redshift (Yoshida et al. 2003). 
For definiteness, we choose M v i r ~ 10 6 M G . To determine the 
pre-burst density, we need to consider the build-up of the H ii re- 
gion around the central Pop III star, and how the density structure 
is modified through the strong photoionization-heating inside 
of it (Whalen et al. 2004). The production of ionizing photons 
strongly depends on the stellar mass, which in turn is determined 
by how the accretion flow onto the growing protostar proceeds 
under the influence of this radiation field (e.g., McKee & Tan 
2008; Hosokawa et al. 2011; Smith et al. 2011; Stacy et al. 
2012). Thus, the assembly of Pop III stars and the development 
of an H ii region around them proceed simultaneously and affect 
each other. The shallow potential wells in the host minihalos, 
with corresponding circular velocities of a few km s _1 , are un- 
able to retain photoionized gas, so that the gas is effectively 
blown out of the minihalo. The resulting photoevaporation has 
been studied with three-dimensional radiative transfer calcula- 
tions (Alvarez et al. 2006; Abel et al. 2007; Greif et al. 2009), 
where one massive Pop III star at the center of the minihalo acts 
as an embedded point source, that also take into account the 
hydrodynamic response of the photoheated gas. 

It is possible to understand the key physics of the photoe- 
vaporation from minihalos with the self-similar solution for a 
champagne flow (Shu et al. 2002). Assuming a p oc r - 2 density 
profile, which describes the typical situation in minihalos out- 
side a nearly flat inner core, one can reformulate the spherically 
symmetric continuity and Euler equations for isothermal gas as 
follows: 
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where x = r/c s t is the similarity variable, and p(r, t) = 
a{x)/4itGt 2 = mnn(r)/X and u(r, t ) = c s v(x) introduce the 
reduced density and velocity, respectively. c s is the sound speed 
of the ionized gas and X = 0.75 the hydrogen mass fraction. The 
density globally decreases with time, after the ionizing source 
has turned on. To derive the immediate pre-explosion value, we 
will therefore set t = f* ~ 3 x 10 6 yr, the typical lifetime of a 
massive Pop 111 star (e.g., Bromm et al. 2001b). 

Different solutions are obtained depending on the ratio 
e = ( c Si i/c s ) 2 , where c s j and c s are the initial and ionized 
isothermal sound speeds. We take T = 3 x 10 4 K and 7) = 7’ vlr , 
where the virial temperature is defined as 


GM m m H 4 / M \ 2/3 /l_ + z\ 

2R wlI k B WmJ V 10 y 


( 3 ) 


The intuition behind the latter identification is that c 2 ; ~ 
k B T vil /m n ~ v 2 k , the specific gravitational potential energy 
of the gas inside a halo, is a measure of how strongly the fluid is 
bottled up, before a champagne flow can occur. Using Table 1 in 
Shu et al. (2002), we can obtain the density profiles at different 
redshifts as follow. First, the value of e can be calculated at 
the given redshift. We then derive the corresponding reduced 
central density, a o, via interpolation. For the given a lt we can 
numerically solve Equations (1) and (2) subject to the inner 
boundary conditions: 


a — ag and v — 0, at x = 0. (4) 



Figure 1. Minihalo circumburst environment. Shown is the hydrogen number 
density vs. distance from the central Pop III star at the moment of its death, 
3 Myr after it was bom. The density profiles are given by the Shu solution 
for different redshifts, as labeled. Typical circumburst densities are ~1 cm~ 3 , 
with smaller values at lower redshifts. The redshift trend reflects the shallower 
minihalo potential wells at later times. 


This relation can be understood as follows. Prior to the 
onset of protostellar collapse, the baryonic density profile in 
the minihalo can be described by the Lane-Emden equation for 
an isothermal sphere 


To avoid the singularity near x = 0, Shu et al. (2002) provide 
the series expansion 


1 d / 2 d(f>\ 
r 2 dr \ dr ) 


= 47TGp 0 e-* la \ 


( 8 ) 


Uq ( 2 \ 2 

a = a 0 + — I - — a 0 J x + ■ ■ ■ , (5) 

2 1/2 \ , 

V = 3 * + 45 \3 _ ° / X + " (6) 

Once the Pop III star has turned on at t = 0, the photoheating 
commences, putting in place a pressurized inner bubble which 
leads to an outward- flowing shock at . The shock obeys 
the usual isothermal jump conditions, written in terms of the 
reduced variables as 

a d 2 

(Vu Xsh)(v d -^sh) — 1? — (^w -^sh) » (7) 

Oiu 


where p B is the density at which the gravitational potential <p 
is set to zero. A simple analytical solution to this equation 
is the well-known singular isothermal sphere (SIS): psis = 
a 2 l(2nGr 2 ), where a 2 — ksTvir/mu. This law approximately 
describes the profile found in realistic simulations of first star 
formation (e.g., Yoshida et al. 2003), outside an inner core 
of radius rj ~ 1 pc, the Jeans length of the primordial gas 
(Bromm et al. 2002). The baryonic mass inside the shock radius, 
r sh = c s t * ~ 100 pc, is about 

rr s h 

Mi = / 47rr 2 psis(z, r)dr. (9) 

Jrj 


where the subscript u(d) indicates the value of the reduced 
velocity and density upstream (downstream) of the shock. 

In Figure 1 , we show the resulting density profiles at the end 
of the Pop III progenitor’s life, and therefore at the completion 
of the photoionization-heating as well. As can be seen, the 
circumstellar densities have dropped significantly from the high 
values, >10 4 cm -3 , present prior to Pop III star formation, and 
are nearly uniform at small radii. Such a flat density profile 
is markedly different from that created by stellar winds. The 
absence of any wind around a Pop III star is indeed predicted, 
as a consequence of very low metallicity (Baraffe et al. 2001; 
Kudritzki 2002), and could serve as an indicator of a low-Z GRB 
progenitor. The post-photoheating density evolves with redshift, 
approximately according to n p ; oc (1 + -), normalized such that 
n ~ 1 cm~ 3 at z ~ 20. This law differs from earlier estimates 
of the Pop III circumburst evolution, where no z-dependence 
was assumed, or a scaling with the density of the background 
universe, n v \ oc (1 + z) 3 (e.g., Gou et al. 2004). 


At the end of photoionization-heating, the mass in the flat inner 
core is about 

47T , 4 tt , 

M f = — r sh p p i - — r^n pi m H X. (10) 

If we neglect any outflows at large radii, we have approximate 
mass conservation, such that M, ~ Mf, resulting in /; pl oc 7’ vlr oc 
(1 + z), as before. 

Our post-photoheating densities are consistent with simula- 
tion results (Kitayama et al. 2004; Alvarez et al. 2006; Abel et al. 
2007; Greif et al. 2009, 2010). We also agree with Whalen et al. 
(2004), if we take into account that their halo had a somewhat 
lower mass, A7 vn ~ 5 x 10 5 M 0 . Empirical support is provided 
in Chandra et al. (2010) who found that the circumburst density 
of GRB 090423 at z = 8.2 was about n ~ 0.9 cm~ 3 by fit- 
ting the radio, X-ray and infrared afterglow, although this burst 
almost certainly did not originate in a minihalo environment. 


3 


The Astrophysical Journal, 760:27 (13pp), 2012 November 20 


Wang et al. 



Figure 2. Atomic-cooling-halo circumburst environment. The density profiles 
are again calculated from the Shu solution, evaluated at t = 3 x 10 6 yr after 
source turn-on. The collapse redshift is here indirectly given by the virial 
temperature, as labeled. We show the case of photoheating from only a single 
Pop III star (solid lines), and that from a stellar cluster (dotted lines). The 
resulting densities again exhibit a nearly uniform inner core, but overall values 
are much higher than in the minihalo case (see Figure 1). 

2.2. Atomic Cooling Halo Case 

Star formation and radiative feedback inside atomic cooling 
halos is significantly less well-understood than the minihalo 
case (e.g., Johnson et al. 2009; Safranek-Shrader et al. 2010, 
2012), and there are no comparable high-resolution simulations 
yet. What is the character of star formation, in terms of IMF 
and stellar multiplicity (Clark et al. 2011a)? Given current 
uncertainty, we here again use the formalism of the Shu solution 
as developed above. We further assume for simplicity that either 
a single Pop III star forms, or a small stellar cluster, with 
properties that would again heat the surrounding Hu region 
to T ~ 3 x 10 4 K. We stress that the numbers derived here 
only provide us with rough guesses at best, and the hope is that 
simulations will eventually become available to firm them up. 

If a cluster of Pop III stars forms, the time where photoheating 
stops will be prolonged. We thus need to evaluate p(r. t) = 
a(x)/4jtGt 2 at 1 = Cluster — U + fsF, where fsF is the timescale 
over which star formation is ongoing, and f* again is the lifetime 
of a single Pop III star. We can estimate ~ f*, because 
disruptive feedback effects will tend to terminate star formation 
when the first member stars die. The resulting densities at the end 
of photoheating are shown in Figure 2. Similar to the minihalo 
case, densities are nearly constant at small radii, but overall 
values are much higher. The reason is that here the potential 
wells are deeper, so that photoheated gas can more easily be 
retained. Typical circumburst densities are n p j ~ 100 cm” 3 . 
Because such large densities would correspond to high afterglow 
fluxes. Pop III GRBs originating in atomic cooling halos may 
be extremely bright, rendering them visible out to very high 
redshifts. For some of the most distant bursts known, such large 
circumburst densities have indeed been inferred. For example, 
Gou et al. (2007) have argued that fitting the GRB 050904 
afterglow at z = 6.3 requires a circumburst density of a few 
100s of cm” 3 . Similar values have been suggested for GRB 
080913 at z = 6.7 (Greiner et al. 2009; Zhang et al. 2009). 

3. MODELING PRE-GALACTIC ENRICHMENT 

In order to predict afterglow spectra from Pop III GRBs, 
we consider explosions that are embedded in a realistic 


cosmological setting. Specifically, we use the first galaxy simu- 
lation carried out by Greif et al. (2010), where the assembly of 
an atomic cooling halo was tracked, resolving all prior Pop III 
star formation in the progenitor minihalos, together with the 
radiative feedback from H ii regions around those stars. In addi- 
tion, the simulation allowed one of the Pop III progenitor stars 
to explode as an energetic SN, explicitly following the transport 
and mixing of the ejected metals into the IGM. We here sum- 
marize a few key aspects of the simulation, and refer the reader 
to Greif et al. (2010) for further details. 

The simulation is carried out in a cosmological box of size 

I Mpc (comoving), and is initialized at z = 99 according to 
a concordance A cold dark matter (ACDM) model with matter 
density Q m = 1 — G\ = 0.3, baryon density Ob = 0.04, 
Hubble parameter h = Hq/(100 km s” 1 Mpc” 1 ) = 0.7, 
spectral index n s = 1.0, and normalization erg = 0.9 (Spergel 
et al. 2003). Density and velocity perturbations are imprinted 
at recombination with a Gaussian distribution, and they are 
propagated to z = 99, when the simulation is started, by 
applying the Zeldovich approximation. Within the highest 
resolution region of our nested, zoomed-in initial conditions, the 
DM and gas particle masses are 33 M Q and 5 Mq, respectively, 
The baryonic mass resolution, roughly given by the mass 
contained inside a smoothed particle hydrodynamics kernel, is 
~400 Mq, close to the relevant Jeans mass in the primordial 
star-forming gas. To capture the chemical evolution of the gas, 
the simulation follows the abundances of H, H + , H”, H 2 , HJ, 
He, He + , He ++ , and e”, as well as the five deuterium species 
D, D + , D”, HD, and HD + . All relevant cooling mechanisms are 
included, i.e., H and He atomic line cooling, bremsstrahlung, 
inverse Compton scattering, and collisional excitation cooling 
via H 2 and HD (Johnson & Bromm 2006). 

In Figure 3, we show the hydrogen number density and 
metallicity averaged along the line of sight within the central 
~ 1 00 kpc (comoving) at two different output times, from z — 
25, briefly after the first star-forming minihalo has collapsed, 
to z = 16.4, closer to the virialization of the first galaxy. 
The distribution of metals produced by the first SN explosion 
is highly inhomogeneous, and the metallicity can reach up 
to Z ~ 10” 25 Zq, which is already larger than any of the 
values discussed for the critical metallicity, Z cr ; t < 1 0 4 Z G 
(e.g., Bromm et al. 2001a; Bromm & Loeb 2003; Frebel et al. 
2007), including those predicted for a dust-driven Pop III/Pop 

II transition, where 1 O' 6 Z Q < Z cr ; t < 1 0 3 Z Q (Omukai et al. 
2005 ; Schneider et al. 2006). The implication is that both Pop III 
and Pop I/II stars will form during the assembly of the first 
galaxies (Johnson et al. 2008; Maio et al. 2010), in turn giving 
rise to the simultaneous occurrence of Pop III and normal GRBs 
at a given redshift (Bromm & Loeb 2006; de Souza et al. 2011). 
Here, we will focus on the case of a Pop III burst exploding in 
one of the (still metal-free) first galaxy progenitor minihalos at 
z — 16.4 (see Section 5). In Figure 4, we show the temperature 
distribution in the vicinity of the explosion site, allowing us to 
calculate (thermal) line broadening in a realistic fashion. 

4. POP III AFTERGLOW EMISSION 
4.1. Basic Physics 

Our modeling of the GRB broadband afterglow emission fol- 
lows standard prescriptions, and we here only briefly summarize 
the basic physics involved (for general reviews, see Meszaros 
2006; Zhang 2007). We specifically consider a relativistic shell 
ejected from the progenitor of a GRB (e.g., Kobayashi 2000). 
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z = 24.68 t H = 130.8 Myr log n H [cm 3 ] z = 24.68 t H = 130.8 Myr log Z [Z e ] 



Figure 3. Possible explosion sites for high-z GRBs. Shown are the hydrogen number density and metallicity contours during the assembly of a first galaxy, averaged 
along the line of sight within the central ~100 kpc (comoving), at two different redshifts. Panel (a): z — 25, briefly after the first Pop III SN exploded. At this time, 
most of the metals reside in the IGM. Panel (b): z — 16.4, closer to the virialization of the atomic cooling halo. Now, metals are being reassembled into the growing 
potential well of the first galaxy. The topology of metal enrichment is highly inhomogeneous, with pockets of highly enriched material embedded in regions with a 
largely primordial composition. These plots are derived from the simulation described in Greif et al. (2010). 


The shell has a rest mass of Mo, an observed thickness of Ao, 
and exhibits an initial Lorentz factor of To = E lso / Mqc 1 , where 
£ iso is the isotropic-equivalent explosion energy. The interaction 
between the shell and the surrounding medium with a density 
of hq produces a forward and a reverse shock. The shocks, in 
turn, accelerate electrons to high energies, distributed accord- 
ing to a power-law: dN(y e )/dy e oc y ( 7 p for y e ^ y m , where 
y e is the internal Lorentz factor for the electrons. The mini- 
mum Lorentz factor is set by assuming that a fraction € e of 
the post-shock energy is transferred to the electrons, resulting 
in: y m = e e y(m p /m e )(p - 2 )/(p - 1), where y ~ T 0 is the 
Lorentz factor of the shocked fluid (e.g., Blandford & McKee 


1976; Sari et al. 1998). If one also assumes that a fraction e# 
of the post-shock energy is deposited into random magnetic 
fields, the gyration of the shock-accelerated electrons around 
the magnetic field will generate synchrotron radiation, thus giv- 
ing rise to the afterglow emission. There is then a second char- 
acteristic Lorentz factor, y c , determining the threshold below 
which electrons do not lose a significant fraction of their en- 
ergy to radiation. Equating the synchrotron cooling time with 
the expansion time of the shell, one finds: y c m e c 2 = P(y c )t ex p , 
where P(y e ) = (4/3)crrcy e 2 (fi 2 /8^) is the synchrotron radia- 
tion power (Rybicki & Lightman 1979), and t exp the expansion 
time, measured in the comoving fluid frame (Sari et al. 1998; 
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Figure 4. Environment for high-z GRBs. Shown are temperature contours 
along the line of sight within the central ~100 kpc (comoving) at z = 16.4, 
corresponding to panel (b) in Figure 3. The center of the growing galaxy is highly 
turbulent, with pockets of cold gas embedded in a hotter, diffuse medium. 


Panaitescu & Kumar 2000). Therefore, the cooling Lorentz fac- 
tor can be written y c = 67zm e c/{crTB 2 t ex p). 

For the resulting spectrum, there are now two cases: in the 
fast cooling case ( y m > y c ), the flux can be calculated from 


F v = iy n 


'/ 3 / v ' 2 


1/3 
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whereas in the opposite, slow-cooling (y m < y c ), case, one has 
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1/3 
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-(P-D/2 
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The break frequencies correspond to the electron Lorentz 
factors introduced above. The frequency v a describes the on- 
set of synchrotron self-absorption, which is important at low 
(radio) frequencies. The afterglow emission exhibits a dif- 
ferent character before and after the reverse shock (RS) has 
crossed the shell. This happens at a time /® ~ (1 + z)Ao/2c = 
16.7 s (1 + ')Ao/ 1 0 12 cm, as measured in the observer frame. 
Before RS crossing (t < f ffi ), we have for the forward 


shock (“/’): 


V/ oc (1 +zr 1 e e 2 4 /2 Ao 3/2 £ i 1 s / 0 2 r 1 , 
v cJ oc (1 +zr 1 e~ V2 nQ 1 A~ l,2 Ery 2 r\ 
v aJ oc (1 +zr 1 4 /S,, o 1/1 ° A o : ll2 E^°r 2 , 

^v.max,/ OC (1 + z)e 1 B /2 n 1 0 /2 E^d^t , 


(ID 


where cIl is the luminosity distance. Here and in the following 
we equate the pre-shock density with that established by 
photoionization-heating, at the time of the Pop III star’s death: 
«o — «pi (see Section 2). For the reverse shock (“r”), we have 

v m ,r OC (1 +z) e e e B r 0 n 0 , 
v c ,r OC (1 +z)- 1 e" 3/2 »0 lA 0 1/2£ i7o 1/2f " 1 - 
<=< (1 4. 

oc (1 + (12) 

After RS crossing, we similarly have 


fm , / — ^(./(l©) J 
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-3/2 
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Vc,f = V c ,/(f®) ( — 

V Iffi 


V a j = V aJ (t®) 


Fv,maxJ — F v ^ m axj(t($) | ^ J , 


for the forward shock. For the reverse shock, we have 

^ —3/2 
-3/2 

J 

- 1/2 


'Vi . r — V mr ft® ) | 

. tea 


v c,r — fc,r(^ffl) I 


V a ,r = V fl , r (r®) ( — 

. tea 


Fy,max,r — F v m ax,rite^) | 

. tea 


(13) 


(14) 


The exact expressions, and further details, can be found in 
Kobayashi (2000) and Wu et al. (2003). Note that the latter 
study assumes the presence of a circumburst wind profile, 
whereas we here consider the constant-density (ISM) case, 
similar to Kobayashi (2000). The basic parameters of the 
afterglow emission region are: To = 300, E IS 0 = 10 53 erg, 
Ao = 10 12 cm, e e = 0.3, = 0.1, and p — 2.5. We have 

verified that our afterglow predictions are consistent with those 
in Gou et al. (2004), and Inoue et al. (2007), as these papers 
focus on Pop III bursts at high redshifts as well. As an example, 
in Figure 5, we show the M- band light curve, corresponding to 
v = 6.3xl0 13 Hz. It can be seen that the reverse shock emission 
dominates during the first few hours, whereas the forward shock 
component takes over at later times. 
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Figure 5. Flux density at v = 6.3 x 10 13 Hz (M band) as a function of observed 
time. We show the emission from the forward shock (dashed line), the reverse 
shock (solid red line), and their combination (solid black line). We use the 
parameters given in the text, for a GRB exploding in a minihalo with circumburst 
density no = 0.5 cm -3 at z = 10. It can be seen that the reverse shock dominates 
during the first few hours, reaching a brightness in excess of ~1 mJy. 
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Figure 7. Observed radio flux (in mJy) at v = 5 GHz as a function of redshift, 
again for a minihalo burst environment. We employ the same convention for the 
lines and symbols as in Figure 6, but we now mark the corresponding sensitivity 
(horizontal line) for the EVLA. As is evident, the radio flux increases with 
observed time, and is well within reach of detection. 
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Figure 6. Observed flux at v = 1.36 x 10 14 Hz (K band) as a function of 
redshift. The lines indicate the flux at different observed times, as labeled. Here, 
we have assumed that the Pop III GRB is triggered in a minihalo environment, 
as modeled in Section 2.1. We mark the A'-band sensitivity of the NIRSpec 
instrument on board the JWST with a horizontal line. The sharp cutoff at z — IV 
is due to Lyot absorption in the IGM. 


4.2. Redshift Dependence 

We here briefly summarize the key aspects of afterglows 
triggered by Pop III stars. Figure 6 shows the observed flux 
at v = 1.36 x 10 14 Hz as a function of redshift for Pop III 
GRBs in a minihalo. According to Section 2.1, circumburst 
densities can then be written as n ~ 1 cm _3 (l +z)/20. The lines 
with filled dots, black triangles and open dots correspond to an 
observed time of 6 minutes, 1 hr and 1 day respectively. The 
straight line marks the A'-band sensitivity for the near-infrared 
spectrograph (NIRSpec) on JWST, estimated for a resolution 
of R = 1000, a signal-to-noise ratio of 10, and an integration 
time of 1 hr (Gardner et al. 2006). Once the observed frequency 
corresponds to an emitted frequency above the Lya resonance, 
v a = 2.47 x 10 15 Hz, which occurs at z ~ 17 for the A'-band, 


all flux will be completely absorbed by the intervening, still 
largely neutral IGM. As can be seen, in the K band the JWST 
will be able to detect GRBs, and to conduct spectroscopy on 
their afterglows, out to z ~ 16 even after 1 day. In the M band, 
the redshift horizon is extended further still, to z ~ 35. At these 
frequencies, the afterglow flux is only weakly dependent on 
redshift. The main reason is that cosmic time dilation implies 
that a given observed time after the trigger corresponds to 
successively earlier emission times, where intrinsic afterglow 
luminosities are much brighter (Ciardi & Loeb 2000; Bromm & 
Loeb 2012). Secondly, circumburst densities modestly increase 
with redshift (see Section 2). 

In Figure 7, we show the flux at v = 5 GHz, again after 
t — 6 minutes, 1 hr, and 1 day. We calculate the sensitivity of 
the VLA, using the following expression (Ioka & Meszaros 
2005): 


F sen 

V 


S/N ■ 2fc B 7sys 

bleffV^tintAv 


23 /rJy 



Ant \~ 1/2 / AP V 1/2 
1 day / \ 50 MHz / 


(15) 


where t, M is the integration time. An the bandwidth, and 
A e ff/T sys ~ 2 x 10 6 cm 2 K 1 . The bandwidth of the Expanded 
VLA (EVLA) 8 can be up to 8 GHz. For a 2 hr integration in 
the 8 GHz band, the EVLA can thus reach a sensitivity of up 
to 6.4 /My, which is shown as a horizontal line in Figure 7. We 
can see that the radio afterglows are easily within reach for the 
EVLA, up to z ~ 40 (e.g., Chandra et al. 2010). 

4.2.1. Peak Fluxes and Wavelengths 

According to hierarchical cosmic structure formation, the 
characteristic halo mass varies as a function of redshift. Thus, 
the typical observed peak flux and wavelength of Pop III GRBs 
are also a function of redshift. In the following, we derive 
ballpark estimates for these key quantities by considering 2 a 
overdensities (e.g., Figure 4 in Clarke & Bromm 2003). This 


http://www.aoc.nrao.edu/evla/ 


7 


The Astrophysical Journal, 760:27 (13pp), 2012 November 20 


Wang et al. 



Figure 8. Observed maximum flux density as a function of redshift at different 
times, t = 1 day, 1 hr and 6 minutes from bottom to top. We assume that the 
Pop III GRB has exploded in a minihalo. As can be seen, afterglow fluxes can 
reach substantial levels early on, and still remain at ~ mJy a day after the trigger, 
out to very high redshifts. 

provides us with the typical halo mass, M xn , at a given redshift. 
In the 2 a case, halo masses become less than 10 s M 0 at 
redshift z > 15. Since this has dropped below the threshold 
mass for Hi cooling (e.g., Yoshida et al. 2003), Pop III stars 
could not form in such low-mass halos. For higher-sigma 
peaks, however, super-critical halos still exist at earlier times. 
Therefore, for simplicity we adopt M wn = 10 5 M 0 at z > 15. 
After the characteristic halo mass is thus determined, the virial 
temperature is fixed as well, and the ratio e = (c s l /c v ) 2 is 
known. The Shu solution then gives the circumburst densities 
at a given redshift (see Section 2), and GRB afterglow spectra 
can finally be calculated. Figure 8 shows the observed peak flux 
at observed time t = 6 minutes, 1 hr and 1 day at z < 30. The 
peak flux varies from millijansky to subjansky. Figure 9 gives 
the observational peak wavelength at the same observed times 
for the same redshift range. The observed peak wavelengths vary 
from millimeter to mid-infrared bands. Early afterglows are best 
observed with facilities such as the JWST, or possibly adaptive- 
optics supported next-generation, extremely large telescopes on 
the ground. At late times, radio facilities provide the optimal 
chance for detection. Among them is ALMA, which will allow to 
probe the afterglow in the millimeter and submillimeter regime, 
where the peak of the emission is located about 1 day after the 
explosion for z ~ 10 sources. 

5. METAL ABSORPTION LINES 

As we have seen above. Pop III GRBs can in principle be 
detected out to very high redshifts. Their afterglow spectra 
provide us with a unique probe into the state of the early IGM, 
including its degree of ionization (e.g., Totani et al. 2006) and 
metal enrichment. The latter will be encoded in absorption lines 
imprinted on the smooth afterglow emission with a signature 
that depends on the nature of the Pop III SNe. Here, we will 
consider for simplicity that prior to the GRB only one nearby 
SN exploded beforehand, dispersing its complement of heavy 
elements into the otherwise pristine IGM. Such a scenario is 
consistent with the clustered nature of Pop III star formation, 
as is expected for the high-er peak host systems (e.g., Greif 
& Bromm 2006). Specifically, we here explore nucleosynthetic 
yields for Type II core-collapse SNe, and for pair-instability 



Figure 9. Observed maximum wavelength as a function of redshift at different 
times, as labeled, for a GRB exploding in a minihalo. The emission maxima 
range from millimeter to mid-IR wavelengths. 


Table 1 

Supernova Yields and Important Transitions 


Element Yx (Scalo) a 

Y x (VMS) b 

Ionization State A., c 

(A) 

f ■ d 

JOSC,l 

c 

0.1 Mq 

4.1 Mq 

Cn 

1334.5 

0.1278 




Civ 

1548.2 

0.1908 




Civ 

1550.8 

0.09522 

0 

0.5 

44 

Oi 

1302.2 

0.04887 

Si 

0.06 

16 

Sin 

1304.4 

0.094 




Si iv 

1393.8 

0.514 




Si iv 

1402.8 

0.2553 

Fe 

0.07 

6.4 

Fen 

1608.5 

0.058 




Fen 

2344.2 

0.114 




Fen 

2382.8 

0.300 


Notes. 

a Metal yield in Scalo IMF case. 
b Metal yield in VMS IMF case. 
c Rest-frame transition wavelength. 
d Oscillator strength. 

supernovae (PISNe). In Table 1, we summarize the average 
yields per SN, Yx , for select elements that are expected to 
dominate the spectral signature. For Type II SNe, we use the 
yields of Woosley & Weaver (1995). These depend on the 
explosion energy, which introduces uncertainties of ~25%. For 
the PISN scenario, we utilize the yields of Heger & Woosley 
(2002, 2010). Table 1 also lists the properties of several key 
transitions, including their rest- frame wavelength, X,, and their 
oscillator strength, f OS c,i- 

We employ a simple model for the ionization structure 
of the metal-enriched gas. If Ihe hydrogen in this region is 
substantially neutral, metals will reside in states typical of H I 
regions in the Milky Way (Cn, Oi, Sin, and Fen for the 
elements of interest), because photons able to further ionize 
these elements will be absorbed by Hi (Furlanetto & Loeb 
2003). The ionization potential of O i is 13.62 eV, and it remains 
locked in charge exchange equilibrium with H i (Oh 2002). A 
modest extragalactic radiation background at frequencies below 
the Lyman limit could easily maintain these metals in a singly 
ionized state (e.g., Furlanetto & Loeb 2003). We thus assume 
that C, O, Si, and Fe are predominantly in the ionization states 
typical for Galactic H i clouds. 
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5.1. Cosmological Radiative Transfer 

Employing comoving coordinates, the equation of cosmolog- 
ical radiative transfer is (e.g., Abel et al. 1999): 


1 dl v ft ■ V/„ 

+ — 

c dt a 


H{t) 

c 



— jv k v I v . 


(16) 


where I v is the intensity of the radiation field, ft a unit vector 
along the direction of the ray, H(t) = a /a the Hubble constant, 9 
and d = ail + dt)/a(t) the ratio of cosmic scale factors separated 
by dt. Here j v and k v denote the emission and absorption 
coefficient, respectively. Equation (16) can be simplified if the 
effects of cosmic expansion are negligible. This will be the 
case if the light-crossing time, hi g ht = L/c, over the region 
of interest, here the (physical) size of our simulation box, 
L = lMpc/(l + z), is small compared to the Hubble time, 
?h — 5 x 10 8 yr [(1 + z)/10]~ 3/2 . Since t\\ g hi hi, the transfer 
equation reverts to its non-cosmological form: 

1 BI V „ 

-~z— +n ■ V/ v = j v — k v I v . (17) 

c dt 

We further neglect any diffuse (re-)emission of radiation, such 
that j v = 0, and assume near steady- state conditions (d/‘dt ~ 0), 
resulting in the static transfer equation: 

h ■ V/ v = —k v I v . (18) 


If 7 U (0) is the intensity of the central afterglow, approximately 
treated as a point source, we have the standard solution: 
I v = /,,(()) exp(— t v ), where r„ is the optical depth along the 
line of sight (LOS). 

We can obtain the optical depth as follows. For a (physical) 
differential distance element, dl, along the given LOS, we have: 
dx v — cr v nxdl, where 

a v = *j7te 2 /(m e c)f OSCii H{u, x)/Ap d (19) 


is the frequency dependent absorption cross-section, and nx the 
number density of species X (e.g., Oi, Cn, Fen, Sin, Hi). In 
the above equation, Avd = bvjc is the Doppler width 10 with 
parameter b = *j2k ti T / The Voigt function, H, is given by 


H(u, x ) 



(x — y) 2 + « 2 


dy. 


( 20 ) 


Here, the relative strength of natural (damping) to Doppler 
broadening is described by the parameter 


r, 

u = , 

4jtApd 


( 21 ) 


where T, is the damping constant for transition i (Morton & 
Smith 1973). The variables x = (v — u,)/Aud and y = v/b are 
the frequency difference relative to the line center in Doppler 
units, and the normalized particle velocity, with respect to the 
GRB location, respectively. The total LOS value of the optical 
depth is then: r v = f cr v nxdl. We note that we can neglect any 
wavelength differences between the GRB source and the metal 


9 We here employ the same cosmological parameters as in the underlying 
simulation (see Section 3). 

10 We here neglect any broadening from turbulent motions. This assumption is 
good in the minihalo progenitors, but breaks down in the emerging atomic 
cooling halo. 



Figure 10. Total GRB spectrum observed at the reverse shock crossing time: 
f® = 16.7 x (1 + 16.4) s. We show the spectral region redshifted into the near- 
IR, accessible to the NIRCam and NIRSpec instruments on board the JWST. 
Metal absorption lines are imprinted according to the Pop III SN event, PISN vs. 
core-collapse (Type II). The former originates from a very massive star (VMS) 
progenitor, whereas the latter from a less massive one. In each case, the cutoff 
at short wavelengths is due to Lya scattering in the neutral IGM. 


absorber, because of the very small expansion-generated redshift 
on the scale of our simulation box. Specifically, we estimate that 
Az = Af/(7.5 x 10 7 yr) ~ 0.005, where we assume z ~ 10, and 
A t ~ hight ~ 3 x 10 5 yr. 

We extract the optical depths along random lines of sight 
with the ray-tracing algorithm of Greif et al. (2009). We 
choose Ng = 160 and N,p = 320 in the spherical coordinates, 
which corresponds to 51200 rays. For each ray, we use 200 
logarithmically spaced radial bins, and approximate the optical 
depth integral with a sum along the LOS. The Pop III GRB 
will likely explode in a still substantially neutral IGM. Any 
flux shortward of the Lya resonance is therefore completely 
absorbed, and because of the corresponding very large optical 
depth, there will be absorption in the red damping wing as well. 
We model this with the analytical formula given by Miralda- 
Escude (1998), using z re i 0 n = 7 and standard cosmological 
parameters (Komatsu et al. 2011). The line shape is not very 
sensitive to the choice of z re ion (Bromm et al. 2001b). 

5.2. Synthetic Spectra 

In Figure 10, we show the part of the resulting afterglow 
spectrum that has been redshifted into the near-IR, where the 
JWST is most sensitive. For the two assumptions regarding 
the Pop III IMF, top-heavy (VMS) and more normal (Scalo), 
the spectra are depicted at the reverse shock crossing time. The 
metal absorption lines are imprinted on the smooth underlying 
GRB afterglow spectrum, rendering their identification straight- 
forward. The cutoff at short wavelengths is due to Lya scattering 
in the IGM which is expected to be still completely neutral at 
z > 10. For detailed predictions, it is advantageous to quan- 
tify the strength of absorption lines through equivalent widths 
(EWs). The observed EWs can be calculated from 

W = (1 + z) /[l- e~ rW ]d\. (22) 

In Figures 11 and 12, we zoom in on the vicinity of the 
strongest lines, again at the reverse shock crossing time for both 
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Figure 11. Spectral signature of metal absorption lines at the reverse shock crossing time: f® = 16.7 x (1 + 16.4) s. We use the metal yields in the VMS case. The 
observed equivalent widths (EWs) are also shown in the figure. If multiple lines are present, the given EWs belong to the strongest lines in each panel. 
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Figure 12. Same as Figure 11, but displaying metal absorption lines for a normal (Scalo) Pop III IMF. Individual EWs are again given for the strongest lines in each 
panel. 


IMF cases. In addition, we explicitly give the derived EWs, 
which are about a few tens of A for a top-heavy IMF, and an 
order of magnitude smaller for a normal (Scalo) IMF. The EW 
values do not change significantly, if evaluated at 1 day after 
the explosion. Such line strengths are well within reach of the 
NIRSpec instrument on board the JWST, which has a spectral 
resolution up to R = X/AX = 1000, and a line sensitivity of 
10“ 18 erg s _1 cm' 2 (Gardner et al. 2006). In general, the limiting 
equivalent width of an unresolved line at observed wavelength 


lobs, at 5er significance, is (Tumlinson et al. 2002): 




51obs 
R( S/N) ’ 


(23) 


where S/N is the signal-to-noise ratio per resolution element. 
Within our assumptions, therefore, spectroscopy of Pop III GRB 
afterglows will be able to distinguish between different enrich- 
ment events, PISN versus core-collapse. There is, however, a 
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Figure 13. Cumulative distribution of metal line equivalent widths. The curves are calculated by tracing 100 randomly selected sightlines through the simulation box, 
centered on the location of the GRB. Top panel: VMS case. The curves show the main lines, as labeled. Bottom panel: normal IMF case, corresponding to conventional 
core-collapse (Type II) SN explosions. It is evident that the resulting line statistics is quite distinct for the two IMF cases. 


Table 2 

Line Flux Ratios from Minimum to Maximum 


Observed Time 
IMF 

Crossing Time 
VMS 

Crossing Time 
Scalo 

1 Day 
VMS 

1 Day 
Scalo 

Oi/Cii 1334.5 A 

0.078-0.086 

0.918-0.965 

0.079-0.090 

0.091-0.097 

Oi/Fen 2334.2 A 

0.014-0.030 

0.408-0.444 

0.014-0.042 

0.040-0.043 

Oi/Fen 2382.8 A 

0.019-0.031 

0.424-0.458 

0.018-0.045 

0.041-0.045 

Oi/Sin 1304.4 A 

0.045-0.059 

0.807-0.881 

0.045-0.089 

0.081-0.088 

Oi/Fen 1608.5 A 

0.078-0.090 

0.800-0.852 

0.078-0.105 

0.080-0.085 


possible complication. The total enrichment level would be 
rather similar to that from one PISN for the case of 10 Type II 
SNe. Such a situation is not implausible, given the number of 
Pop III star formation sites, i.e., minihalos, in the vicinity of 
the Pop III GRB (see Section 3). However, we could then still 
discriminate between the cases using absorption line flux ratios, 
which remain quite different. 

In Table 2, we list ranges for select flux ratios, as encountered 
along different LOSs. We consider the ratios Oi/X, where X 
denotes another species, at both the reverse shock crossing time 
and at 1 day after the explosion. In the former case, there 
is an order of magnitude difference between the two IMFs, 
allowing a robust identification. At later times, however, this 
clear-cut signature disappears again. In Figure 13, we show the 
cumulative EW distribution for 100 randomly selected sightlines 
at redshift z = 16.4. Again, the resulting distributions are quite 
different for the two IMF cases. It is, however, important to 
note that the statistical significance of this analysis is limited by 
the fact that it is based on only one simulation, and only one 
GRB afterglow. In reality, observations will sample multiple 
explosion sites with respective LOS distributions that sample 
many clouds at various evolutionary stages along the way. Our 
analysis should still provide a representative view of the typical 


situation, whereas a more complete study has to await the further 
increase in available computational power. 

6. SUMMARY AND CONCLUSIONS 

In this paper, we develop a diagnostic to study pre-galactic 
metal enrichment in the vicinity of the first galaxies. Specifically, 
we utilize the bright afterglow of a Pop III GRB as a featureless 
background source, and calculate the strength of metal absorp- 
tion lines that are imprinted by the first heavy elements produced 
by Pop III SNe. To approximately derive the metal absorption 
line statistics, we use an existing highly resolved simulation of 
the formation of a first galaxy which is characterized by the on- 
set of atomic hydrogen cooling in a halo with virial temperature 
>10 4 K. The reverse shock initially dominates the afterglow 
flux a few hours after explosion, followed by the forward shock 
emission later on. The fluxes in the near-IR and radio bands 
may be detectable with the JWST and the VLA out to z > 30. 
We predict that the afterglow emission peaks from near-IR to 
millimeter bands with peak fluxes from mJy to Jy at different 
observed times. Metal absorption lines in the GRB afterglow 
spectrum, giving rise to EWs of a few tens of A, may allow us 
to distinguish whether the first heavy elements were produced 
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in a Pop III star that died as a PISN, or a core-collapse SN. To 
this extent, the spectrum needs to be obtained sufficiently early, 
within the first few hours after the trigger. 

The absorption signature of the first SN events might allow us 
to constrain the underlying mass scale of Pop III stars. This is the 
key input parameter to predict their evolution, nucleosynthetic 
yields, and modes of death, which in turn determine how the 
first stars impact subsequent cosmic history. The strength of 
the associated stellar feedback governs the assembly process of 
the first galaxies, in the sense that the stronger feedback from 
more massive stars shifts their formation to later stages in the 
hierarchical build-up of structure (e.g., Ricotti et al. 2002; Greif 
et al. 2010; Frebel & Bromm 2012; Ritter et al. 2012; Wise et al. 
2012). 

It will, however, be very challenging to directly probe the 
initial epoch of cosmic star formation. The reason is that even the 
JWST will not be able to detect individual first stars, but instead 
only more massive stellar systems or clusters that form later on 
in more massive host systems (e.g., Pawlik et al. 201 1). Pop III 
GRBs may thus afford us one of the few direct windows into 
the crucial epoch of first light, another one being the extremely 
energetic pair-instability SN or hypernova explosions that are 
predicted to end the lives of the most massive Pop III stars (e.g., 
Pan et al. 2012; Hummel et al. 2012). The promise provided by 
the combination of a wide field GRB trigger mission, such as 
JANUS or Lobster, with the next-generation of highly sensitive 
near-IR telescopes, such as the JWST or the planned ground- 
based extremely large facilities (the E-ELT, 11 GMT, 12 and 
TMT 13 ), is huge. GRBs are likely to play a key role in finally 
opening up the high-redshift frontier, all the way back to the 
very beginning of star and black hole formation. 
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